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VACUUM EXPECTATION VALUE ASYMPTOTICS 

FOR SECOND ORDER DIFFERENTIAL 
OPERATORS ON MANIFOLDS WITH BOUNDARY 



Thomas P. Branson^, Peter B. Gilkey*, Dmitri V. Vassilevich* 



Abstract. Let M be a compact Riemannian manifold with smooth boundary. We 
study the vacuum expectation value of an operator Q by studying Tr£2 Qe~*^, where 
D is an operator of Laplace type on M, and where Q is a second order operator 

^ ' with scalar leading symbol; we impose Dirichlet or modified Neumann boundary 

Cu , conditions. 



> 

^yQ . §1 Introduction 

, ' Let M be a compact smooth Riemannian manifold of dimension 771 with smooth 

fvq . boundary dM. We say that a second order operator D on the space of smooth 

f^ ' sections C°°(V) of a smooth vector bundle over M has scalar leading symbol if 

C^^ , the leading symbol is h^-' Iv £,i£,j for some symmetric 2-tensor h. We say that D is 

^1* ' of Laplace type if /i'-' is the metric tensor on the cotangent bundle. Let D be an 

operator of Laplace type. If the boundary of M is non-empty, we impose Dirichlet 

or Neumann boundary conditions B to define the operator Dg, see §4 for further 

Q-i. details. Let Q be an auxiliary second order partial differential operator on V with 

~H ' scalar leading symbol; if the order of Q is at most 1, then this hypothesis is satisfied 

trivially. As t J, 0, there is an asymptotic expansion 



X 

c3 ■ "=-2 



(1.1) Tri2(Qe-*^'')^ ^ an{Q,D,B)t 



(n-m)/2 



see Gilkey [8, Lemma 1.9.1] where a different numbering convention was used. The 
invariants an{Q,D,B) are locally computable. We have a-2{Q,D,B) — and 
a^i{Q,D,B) = if Q has order at most 1. If the boundary of AI is empty, the 
boundary condition B plays no role and we drop it from the notation; in this case, 
if n is odd, then a„(Q, D) — 0. 
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2 T. BRANSON, P. GILKEY, AND D. VASSILEVICH 

Our paper is motivated by several physical examples. First, consider a Euclidean 
quantum field theory with a propagator D~^ depending on external fields. Typ- 
ically, _D is a second order differential operator of Laplace type. In the one-loop 
approximation, the vacuum expectation value < Q > of a second order differential 
operator is given by < Q >= Tr]^2{QD^^). By formal manipulations, this can be 
represented in the form 



<Q>^ S;^dt TrL2{Qe-'D) ^ ^ dtj:Z-2^niQ,D,B)t(^- 



m)/2 



These integrals are divergent at the lower limit and need to be regularized. This can 
be done by replacing by 1/A in the limits of integration; A is called the ultraviolet 
cutoff parameter. The coefficients an{Q,D,B) define the asymptotics of < Q > as 
A — > cxD. The first m terms are divergent and are essential for renormalization. The 
coefficients an{Q, D, B) also define large mass asymptotics of < Q > in the theory 
of a massive quantum field; for details see for example [1]. 

A second example is provided by quantum anomalies. In the Fujikawa ap- 
proach [6], the anomaly A is defined as ^ = \iin\^oo Tr{Qe~^^^ ), where Q is 
the generator of an anomalous symmetry transformation, and D is a regulator. 
Usually divergent terms may be absorbed in renormalization and one has that 
A ^ am-2{Q,D,B). Other examples where these asymptotics arise naturally are 
the study of the anomaly for an arbitrary local symmetry transformation, and in 
the study of the vacuum expectation value of the stress-energy tensor. 

In this paper, we will study the asymptotics a„(Q, D, B) in a general mathemat- 
ical framework. In §2, we review the geometry of operators of Laplace type and 
derive some variational formulas. The operator D determines the metric 5, a con- 
nection V on V , and an endomorphism E of V . Conversely, given these data, we 
can define an operator of Laplace type D{g, V, E)\ see Lemma 2.3 for details. Let 
q2 be a symmetric 2-tensor and let qi be a 1-form valued endomorphism of V . We 
shall use 52 to define a variation of the metric g{e) :~ g + eq2 and we shall use qi to 
define a variation of the connection V(e) := W + eqi. Let Q2 :— deD{g{e),\7, E)\^=o 
and let Qi :— dsD{g,V{e), E)\g=o. Let Q be a second order operator with scalar 
leading symbol. We may decompose Q = Q2 + Qi + Qo for Qq G End(V^) and for 
Q2 and Qi defined by suitably chosen q2 and qi. Since an{Q,D) — '^^an{Qi,D), 
it suffices to compute the an{Qi, D). In Lemma 2.4, wc will study the operators 
Q2 and Qi and show that 9ga„+2(l, -D(g),S) = —an{dgD{g),D{Q),B) for any 1- 
parameter family of operators of Laplace type and fixed boundary condition B. In 
§3 and §4 we use this variational formula and apply results of [2] and [5] to study 
the invariants an{Q, D,B); in §3 we consider manifolds without boundary and in 
§4 we consider manifolds with boundary. 

An operator A is said to be of Dirac type if A^ is of Laplace type. Branson and 
Gilkey [2] studied the asymptotics of Ttl2 {Ae^*"^ ) for an operator A of Dirac type 
on a closed manifold. In §5, we use the results of §3 to rederive these results and 
to compute some additional terms in the asymptotic expansion. The numbering 
convention we shall use in this paper differs from that used in [2]; the invariants 
a„(j4, A^) of this paper were denoted by a„-i(A, A^) in [2]. 

§2 Geometry of operators of Laplace type 

We adopt the following notational conventions. Greek indices /i, v, etc. will 
range from 1 through m = dim(M) and index local coordinate frames 9^ and dx'^ 
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for the tangent and cotangent bundles TM and T*M. Roman indices i, j will also 
range from 1 through m and index local orthonormal frames e^ and e* for TM and 
T*M. We shall suppress the bundle indices for tensors arising from V. We adopt 
the Einstein convention and sum over repeated indices. Let D be an operator of 
Laplace type. This means that we can decompose D locally in the form 

(2.1) D = -{g''^Ivd,d^ + a''d, + b) 

where a and b are local sections of TM (g) End(F) and End(y) respectively. It is 
important to have a more invariant expression than that which is given in equation 
(2.1). Let r be the Christoffel symbols of the Levi-Civita connection of the metric 
g on M, let V be an auxiliary connection on V, and let E G C°° (End{V)) . Define: 

.2 2^ D{g,V,E):=-iTTgV' + E) 

We compare equations (2.1) and (2.2) to prove the following Lemma. 

2.3 Lemma. // D is an operator of Laplace type, then there exists a unique con- 
nection V on V and a unique endomorphism E of V so that D = D{g, V, E). 

(1) If Lu is the connection 1-form ofV, then lus = giysii" + 5^'^r^o-''/y)/2. 

(2) We have E = b- .g'"^(a^w^ + uj^ujf, - t^^L^^'"). 

Let D = D{g, V, E). We use the Levi-Civita connection of the metric g and the 
connection V on V to covariantly differentiate tensors of all types. We shall let ';' 
denote multiple covariant differentiation. Thus, for example, Df = —{f-kk + Ef). 

2.4 Lemma. Let D = D{g, V, E) be an operator of Laplace type. Let q2 = q2,ij be 
a symmetric 2-tensor and let qi = qi^ be an endomorphism valued Ttensor. Then 

(1) Oi/ := d,D{g, V + equE)fl=o - -2gi../;. - ql,^,^f . 

(2) Qif := deD{g + £q2, V, E)f\s=Q = 12.13 f,i] + (2g2,jju- - g2jj;i)/;j/2- 

(3) Let D{q) be a smooth I -parameter family of operators of Laplace type and 
fix B. Then a„{dgD{g), D{g), B) = -dgan+2{l,D{g),B). 

Proof. Fix a point xq G M; we may assume that xq is in the interior of M. Choose 
coordinates centered at xq and a local frame for V so that ^^^(xo) = (5^^, r(a;o) = 0, 
and so that uj{xo) = 0. We use equation (2.2) to compute: 

d,D{g, V + eqi,E){xo)U=o = -g'"'{2qi,^d^ + d,qi^„){xo) 
= (-2gi_iVi -~qi.i-i){xo), 

deig + £q2Y'^ {xo)\e=0 = -q2.Ma{xQ), 

2deT{g + Sq2)„a'^{xo)\s=0 = {(l2,tj.u;a + q2,t^a-u - q2,ua-ii){xo) , 

deD{g + eq2,V,E){xo%=o 

= {q2,iya{duda- + d^UJ^) + deT{g + £q2) ^if," du)(xo%=Q. 

The first two assertions now follow. We use [8, Lemma 1.9.3] to see that the 
asymptotic series of the variation is the variation of the asymptotic series. We 
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equate coefficients in the following two asymptotic expansions to complete the proof: 

-^ - Efc ak{d,D{g), D{g),B)t(^+^-"^)/\ D 

To use Lemma 2.4, we shall need some variational formulas. Let R^ua^ be the 
curvature of the Levi-Civita connection with the sign convention that the Ricci 
tensor is given by pi,^ := Rfj,ua^ and the scalar curvature is given by t := g^'^ Pua- 
Let Aq = 5d be the scalar Laplacian, let cfvol be the Ricmannian measure, and let 
_F^i/ be the curvature tensor of V. 

2.5 Lemma. Let V(e) := V + eqi and g{g) := g + gq2. Let Tij := (de\e=oF)ij, 
Ti^ivJ := {dg\g=oR)iJ^ucy^ , and let V := dg\g=o/!^o- Then 

(1) Tij = gij;i - qis-j and 9e|e=o(VF)ij;fe == Tij-k + [gi,fc,.Fy]. 

(2) Vf = q2,ijf;ij + (2(72, y;j - q2jj:i)f;l/2. 

(3) dgdvo\\g^o ^ q2Aidvol/2. 

(4) (dgT)f,i,''\e=o = g"'' {q2.,t,r-M + 92,1.7;^ - q2,t,u;-,)l'^- 

(5) TZ^iya ~ g ^ \q2,na:iv + q2,vi;an ^ 'Z2,/J7;cri^ ^ '?2,!/cr;7/i ^ 92,cr/3-R/ii^7 

^92,7/3-R^iy(T )/2. 

(6) dg\g=QT = —q2,ijPij + Tlkiik- 

(7) 9g|g=o|p| ~ '^T^kijkPij — 2q2.ijPikPjk- 

(8) Sglg^ol-Rl = ^TZijklRijkl — '2q2,jnRijklRinkl- 

Proof. The assertion (1) is immediate from the definition. Assertion (2) follows 
from Lemma 2.4. Assertions (3) and (4) are straightforward calculations. Assertion 
(5) follows from assertion (4) and from the identity^ : 

Q2,(T^:L'fi Q2,(T^:fi!y — q2.'yp^^u(7 q2.(7p^fj.u^ • 

Raising and lowering indices does not commute with varying the metric so we 
emphasize that the tensor TZ is the variation of a tensor of type (3, 1). The remaining 
assertions now follow. D 

§3 Manifolds without boundary 

Lemma 2.4 reduces the computation of a„(Q, D) to the special cases an{Qi, D) 
for i = 0, 1,2. Recafi that an{Q,D) = forn odd; a-2{Q,D) = if ord((3) < I. 
If P is a scalar invariant, let V[M] := Jj^ P{x)dvol{x). Let try be the fiber trace. 
We refer to Gilkey [7] for the proof of the following result: 

3.1 Theorem. Let M be a compact Ricmannian manifold without boundary, let 
D be an operator of Laplace type, and let Qo G End(y ) . Then 



(1) ao{Qo,D) = {4Tr)-"^/Hrv{Qo}[M]. 



^We are grateful to Arkady Tscytlin who pointed out a sign error in this identity in the previous 
version of the paper 
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(2) a2{Qo,D) = (47r)-'"/26-itrv{Qo(6S + r)}[M]. 

(3) a4(Qo, D) = (47r)-'"/2360-i try{Qo(60£;;fefc + 60rS + ISQE^ 

+ l2T,kk + 5t2 - 2|p|2 + 2|i?|2 + 30F,,F,,)}[M]. 

(4) ae{Qo.D) = [A^)-"^'^ trv{QolV.{l%T,un + nT,kT,k - 2p,,,uPir,k 

+24pjkPj„-kn + 12 RijkiRtjke-nn + 35/9r^ - lA/Srp^ 

+ 14/3Ti?^ - 20S/9pjkPjnPkn - ^-^/'iPijPklRikjl 

— W/SpjkRjnuRknei — 4:4:/9RijknRijipRknep 

—80/9RijknRiekpRjinp) + 360 Qoi^iFij-kFij-k + 2Fij-jFik-k 

'^-^ij'.kk-'^ij I ^^ ij ^ ij-.kk ^■^^ ij ^ jk^ ki ^^ijk7i^ ij ^ kn 

—4:PjkFjnFkn + ^TFknFkn + (iE-ajj + SOEE-a + SOE-uE 
+30E.^E.^ + 60E^ + UEF^jFij + QF,jEF,j + UFi^F^jE 
+ mTE.kk + ■ipjkE.jk + l2T,kE,k - 6E.jFij,i + GFij.iE.j 
+30EET + l2ET,kk + 5Et^ - 2EpjkPjk + 2ER^)}[M]. 

Next, we study the invariants an{Qi, D). 

3.2 Theorem. Let M be a compact Riemannian manifold without boundary, let 
D be an operator of Laplace type, and let Qi = d^D{g,V + eqi, E)\^^q. Then 

(1) ao(Qi,^) = 0. 

(2) a2{Qi,D) = -(47r)-"/2360-itry{60Fy.F.j}[M]. 

(3) ai{QuD) = -(47r)-™/2360-itry{-8i^y;fe^y;fe -8Fy-;fegi,feFy- 

+8Fij-kFijqi^k + 4:Fij-jTik:k + 4:Fij-jqi,kFik — AFij-jFikQi^k 

'J^-^ij ^ jk^ki ^■^^ijkn^ ij-^kn ^Pjk^jn-^kn ^ ^'^^^kii-^kn 

-GOE.Mi.kE + 60E,kEqi,k + 30EF,jT,j + SOFT,, F,j}[M]. 

Proof We use Lemma 2.4 and Theorem 3.1. Note that Ao(/) = ^f-kk is indepen- 
dent of the connection V for / G C°°{M). We compute: 

(1) d,tTv{60E,kk)U=o = 60d,tTv{F);kk\e=o - -QOd,AotTv{F) = 0. 

(2) detTv{30F^)\e=o - trv(60F,,.F„). 

(3) 9e tTv{8Fij.kFij-k + QFij-kkFij + 6FijFij.kk)\e=o 

= 9e{try(-4Fij;feFy;fe) + trvi6FijFij).kk}\e=o 

— tii:v{~'S>Fij-k^ij;k — S>Fij-kqi,kFij + 8Fij-kFijqi,k)- 

(4) d, trv{30EE,,, + 30E,,,E + 30E,,E,,)U=o 

= d,{trv(-30E,,E,,) + tiv(30EE),,,}U=o 
^trv{-Q0E,kqi,kE + 60E,kEqi,k). D 

We conclude this section by studying an{Q2,D). The foUowing resuh is a con- 
sequence of Lemmas 2.4, Lemma 2.5, and Theorem 3.1. We omit the formula for 
0'4{Q2, F)) in the interests of brevity. 
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3.3 Theorem. Let M be a compact Riemannian manifold without boundary, let 
D be an operator of Laplace type on C°^{V), and let Q2 = deD{g + eq2,y , E)\^=o. 

(1) a_2(Q2,-D) = -iflo ((72, ii,-D). 

(2) an{Q2,D) - -^a2{q2.u,D) ^ {^^T)-"'/^6-HTv{{-q2,^JP^J +nkuk)Iv}[M]. 

(3) a2{Q2,D) ^ -\ai{q2M,D) ~ (4^)-™/2360-i{-Ptry(60S + 12T/y) 

+ try{10( — g2,ijPy + T^kiik)Tlv — 2{2TlkijkPij — 2q2.ijPikPjk)Lv 
+2{2TZijklRijkl — 2q2JnRijklRinkl)Iv + (^0{ — q2,ij Pij + T^kiik)E 
-60q2^^kF^JF,k}[M]. 

§4 Manifolds with boundary 

We now suppose M has smooth non-empty boundary dM. Near dM, let Ci be 
a local orthonormal frame for TM where we normalize the choice so that Cm is the 
inward unit normal. We let indices a,b, ... range from 1 through m — 1 and index 
the resulting orthonormal frame for the tangent bundle T{dM) of the boundary. 
Let / e C°°(V). Let S be an endomorphism of V defined on dM. The Neumann 
boundary operator is defined by B^f := {'Vm + S)f\QM and the Dirichlet boundary 
operator is defined by Bg f = f\dM] we set S ~ with Dirichlet boundary condi- 
tions to have a uniform notation. Let Lab ~ {^Ca^b, &m) be the second fundamental 
form. Let ':' denote multiple covariant differentiation tangentially with respect to 
the Levi-Civita connection of the metric on the boundary and the connection V 
on V. The difference between ';' and ':' is given by the second fundamental form. 
For example, E.^a and E.,a agree since there are no tangential indices in E to be 
differentiated while E.^ab = E.,ab — LabE-^m- There are some new features here which 
are not present in the case of manifolds without boundary in the following formulas. 
The invariants an{Q , D , B) are non-zero for odd n and the normal derivatives of Qq 
enter. We still have a_2(Q, D,B)^0 if ord(g) < 1. 

4.1 Theorem. Let M be a compact Riemannian manifold with smooth boundary, 
let D he an operator of Laplace type and let Qo G End(y). Then 

(1) ao{Qo,D,B^) = (4^)-"/2try{Qo}[M]. 

(2) ai(go,i?,B|) = ±(47r)-("-i)/24-itry{Oo}[9M]. 

(3) 02(^0, ^,^1) = i'iTr)-"'/^6-^{trv{Qo{&E + t)}[M] 

+ trv{Qoi2Laa + 125) + (3+, -3-)Qo:m}[9M]}. 

(4) aaCQo, D, Bf) = (4^)-('"-i)/2384-i trv{Oo((96+, -96-)^ + (16+, -16~)t 

+ (8+, ^8-)Ramam + (13+, -7-)LaaLbb + (2+, 10-)LabLab + QQSLaa 

+192S^) + Qo:m((6+, 30-)L„a + 965) + (24+, -24-)Qo-m„^}[dM]. 

(5) ai{Qo,D,B'l) = i4TT)-"'/^360-^{tTv{Qo{Q0E.^kk+&0TE+l80E^+30FijF,j 

+ l2T.,kk + 5t2 - 2p2 + 2i?2)}[M] 

+ try{Oo(240+, -120-)S;„ + (42+, -18-)t;„ + 24L„„:b6 + UOELaa 

+ 20TLaa + 4:RamamLbb — ^2RambmLab + '^RabcbLac + 360(5-© -|- ES) 

+21-^ {{2%Q+ ,m-)LaaLbbLcc + (168+, -2M-)LabLabLcc 

+ (224+, 'i2{r)LabLbcLac} + 1205T + lUSLaaUb + '^^SLabLab 

-h48052L„a + 48053 + 1205aa + Qo;m((180+,-180-)£;+(30+,-30-)r 
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+ (84+, -180-)/7 • LaaLbb + (84+, 60-)/7 • LabLab + 72SLaa + 240^2 
+ Q0;rnmi2iLaa + 1205) + (30+, -30-)0o;«»m } [9M] } . 

4.2 Theorem. Let M be a compact Riemannian manifold with smooth boundary, 
let D be an operator of Laplace type, and let Qi — di;D(g, V + eqi, i?)|e=o- Then 

(1) a_i(Qi,i?,S±) = 0. 

(2) ao{Qi,D,B^) = -(47r)-"/26-i try((-12+, 0-)gi,,„)[9M]. 

(3) ai{Qi,D,B^) = -(384)-i(47r)-("-i)/2try((-96+,0-)gi^„L,a 

^384Sqi,„,)[dM]. 

(4) a2{Qi,D,B^) ^ -(47r)-"/2360-i{try(60Fy.Fy)[M] 

+ trv{{-720+,0-)qi^,nE+{-120+, 0-)gi,™T+(-144+, 0-)gi.„,Laaifc(, 
+ (-48+, 0-)qi.^LabLab - 960Sqi.^Laa - 14405V,m)[aM]}. 

(5) a3{Qi,D,B-) - 5760-i(47r)("-i)/2try{240F,b^,b - 720 FamTam}[dM]. 

(6) // i/ie boundary of M is totally geodesic, then 
a3{Qi,D,Bf) = -5760-H4^)(™-i)/2trv(-1440^;„gi,„ 

+ 1440(gi^,„S - Eqi^rn)S + 240FabJ'ab - 960TSqi,m - 240p™„5qi,„ 

+ 180FamJ^am — 270T-mqi,m + 720 S;aqi,m:a 

+720S..a{qi,aS - Sqi.a) - 2880E{Sqi.m + qi.^mS) - 5760qi.,nS^)[dM] . 

Proof. If Qo = qolv for qo G C^{M) is a scalar operator, then Theorem 4.1 follows 
from Branson and Gilkey [2]. If Qo is not a scalar operator, we must worry about 
the lack of commutativity; the only point at which this enters is in the coefficient 
of tiviQoSE) and tvviQoES). We express 

aiiQo, D, Bs) = (47r)-"/23go-i tvv{CiQoSE + C2QoES)[dM] + other terms; 

the sum Ci + C2 = 720 is determined by the scalar case. If D, Qqj and S are real, 
then Tr2^2 (Qoe~*^) is real; this shows that Ci and C2 are real. If Qo D, and S are 
self-adjoint, Tri2(Qoe-*^) is real so tYv{Qo{CiSE + C2ES))[dM] is real; this now 
shows Ci = C2 and completes the proof of Theorem 4.1. 

To keep boundary conditions constant, let S{e) := S — sqi^m so deS\e=o = —qi.m- 
Assertions (l)-(5) of Theorem 4.2 now follow directly from Lemma 2.4, from Lemma 
2.5, and from Theorem 4.1. In [5], we showed that 

05(1, D, S|) = ±576G-i(47r)('""^^/2 try{(36G£;;„„ + 144G£;;„S' 
+72G£;2 + 240£;:aa + 240t£; + UOFabFab + 48r;,i + 20t2 - 8p^ 
+8i?2 _ l20pmrr.E - 2Gp„„r + 48Gr52 + (9G+,' -360-)F„„F„„ 

-\-bijijO:a^:a + ^^J^ammbPah ~ ^'PmrnPmrn ~ ^^J^ammbJ^arnmb + ZooU-C/O 

+144G5'4) + £:}[aM] 

The variation of the terms other than £ gives rise to the expressions listed in 
Theorem 4.2 (5,6). The remainder term £ is given below. It vanishes if the boundary 
is totally geodesic and involves 40 undetermined coefficients. 

O ^ Cl-^ Ijaa-t^-.ni 'y Q>2 ^aa'^-.m 1 ^3 ^ab-t^animb-.ni ^" ^4 ^aa^:bb i 0(5 J^ab^:ab 
+"6 Eaa:b^:b + dj Lab:a^:b + Sg ^aa:bb^ + "9 Eab-.ab^^ + d-^QLaa-.b^cc-.b 
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~^^ll^ab:a-^cc:b y ^i2^ab:a^bc:c i ^i^^ab:c^ab:c i ^i4-^ab:c-^ac:b 
'T-d-^r^Laa:bb^cc~rdiQhab:ab^cc~rdi'jLab:ac^bc~rdij^Laa:bc^bc~rdig-Lf,caaJ-'bc 
+ UAO+LaaSE + dt^LaaSpmra + 240+LaaS'T + d+iLabSpab 
+dt2LabSRmabm + (195+, lQb-)LaaLbbE + (30+, ISO" )iafciah-E 
+ (195 + /6, IQ^~ /&)LaaLbbT + i5+,25^)LabLabT + d^^LaaLbbPmm 
'^d24J^ab^abPmni r d2!^J^aa^bcPbc i d2^1^aa^bc^mbcm i d2jJ^ab^acPbc 
~\-d2g,J^abJ^acJ^mbcm + d2QJ^ ab^ cd^acbd + d^^Laa^^ + d^^Laa^bb^^ 
-\-d^2^abi-'abi3 + d^^Laa^bb^cc^ + d^^LaaEbcJ-'bc^ + d^^Lab^bcLacJ 
+d^QLaaLbbLccLdd + d^j L aaLbbL cdL cd + d^^L abL abL cdL cd 
+d^gLaaLbcLcdLdb + d^f^LabLhcLcdLda 

The variation of £ is zero for Dirichlet boundary conditions or if the boundary is 
totaUy geodesic. D 

To study an{Q2, D,Bg) we need some additional formulas. 

4.3 Lemma. 

(1) Let g{g) := g + gq2, M := 9ele=oem(e), anrf /Iq/j := 5e|g=oia/3- T'Aen 

a) A/" = -(72, amCa - <?2,mmem/2. 

b) /3ah = {q2,am;b + <l2,bm;a — <l2Mb;rn — q2,mmLab) f^. 

(2) q2,am;a = q2,am:a ^ -^00^2, mm + Labq2,ab- 

(3) T^-fciifc = q2M-M ~ q2,ii;kk- 

Proof. Let 1 < a,/? < to — 1. Let y = (y") be local coordinates on the boundary 
dM centered at j/o- We suppose gapivo) = ^q/3 + 0{\y\'^). Introduce coordinates 
X — (jj,x™) so the curves t h^ {y,t) are unit speed geodesies perpendicular to the 
boundary. Then gmm — 1 and gam — 0; dm is the inward geodesic normal vector 
field for g. Let N{g) be the inward geodesic normal vector field for the metric g{g). 
Expand N{g){yo) = dm + Q{c"^dm + c'^^^) + Olg^). We prove the first assertion by 
solving the equations 

=9{Q){N{g),da){yn) = g{d^ + q2,c.m){yo) + Oig") 

1 =g{Q){N{g).N{g)){yo) - 1 + g(2c™ + g2,mm)(yo) + 0{g^) 

to see c"(j/o) = —q2,ma{yo) and c'^{yo) = —(72, mm/2. We use Lemma 2.5 to compute 
the variation of the Christoffel symbols and complete the proof by computing: 

Lap^Tip)ap'g{N{g),d^) 

^afliVo) =(ra/3™ ~~ q2,mmLa(i/'2){ya)- 

The second assertion is immediate, the third follows from Lemma 2.5. D 

Dirichlet boundary conditions are unchanged by a variation of the metric g. The 
following result follows from Lemma 2.4, Theorem 4.1, and Lemma 4.3. We omit 
the formula for 02 in the interests of brevity. 

4.4 Theorem. Let M be a compact Riemannian manifold with sm,ooth boundary, 
let D be an operator of Laplace type, and let Q2 '■— de\e=f)D{g + £q2^ V, E). Then 

(1) a-2{Q2,D,B-) = -(47r)-"/2tryte,,»/2}[M]. 
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(2) a^i{Q2,D,B-) = (4^)-("-i)/24-itry{<72.aa/2}[aM]. 

(3) ao{Q2,D,B-) = -(47r)-"/26-i{try{g2,n(6S + r)/2 - g2,^JPy 

+7^fe^^fc}[M] +tlv{q2,aaLaa + 2Caa - 2q2,abLab}[dM]} . 

(4) ai(Q2, D, B-) ^ (4^)-("-i)/2384-i try{(?2,aa(96^ + 16t + 8Ra 

+ 7LaaLbb — WLabLab)/'2 + 16{ — q2,ijPij + TZkiik) 

+ 8[R-amam ^ <l2,abRambm + '^RamaAfJ +7{'2^aaLbb ~ '2<l2,abLabL cc) 

-lQ{2LabLab - 2q2,acLabLcb)}[dM]. 

(5) Suppose for simplicity that the metric g is flat, i.e. that Rijki = 0. Then 
a2{Q2,D,B-) = -(47^)-"/2360-l{-X>try(60£;) + try(607^iJJ,£; 

+30g2,«^;fefc + 90(72,mS' + 15q2.un^)}[M] 

-(4^)-'»/2360-ltry{60S;„(72,m™ + 120E,aq2.arn - 187^,fcfc,;„ 

'T'J-^aa\'^^'^ikki y ^''^bmbm) ~r ^'^'^ambm-'-^ab ~r '^'^abcb-'-^ac 

+ q2M-^(^E.,m + GOELaa + l2Laa:bb 

+ 20/21LaaLbbLcc + U/7LabLabLcc + 160 /21LabLbcLac) 

{120E + AO/7 LaaLbb + 88/7 LabLab){Ccc - q2.cdLcd) 

+ 88/7{2£abLabLcc — 2q2.bdLabLadLcc) 

+ 320/7{£abLbcLca — q2,adLabLbcLcd) + 12Lbb;cq2,aa:c}[dM]. 

We now study Neumann boundary conditions. The situation is quite different 
as Neumann boundary conditions are not invariant under general perturbations of 
the metric; if qam 7^ on dM , Bg{g) wiU involve tangential derivatives regardless 
of how 5* is varied. Thus Lemma 2.4 is not directly applicable. Nevertheless, we 
can still compute the first three terms in the asymptotic expansion. 

4.5 Theorem. Let M be a compact Riemannian manifold with smooth boundary, 
let D be an operator of Laplace type, and let Q2 '.— deD{g + eq2, V, E)\^^q. Then 

(1) a^2{Q2,D,B+) = -i4Tr)-"^/Hrv{q2,u/2}[M]. 

(2) a^,{Q2,D,B+) - -(4^)-(™-i)/24-itry{g2.aa/2}[9M]. 

(3) ao{Q2,D,B+) = -(4^)-'"/26-i try{(?2,^^(6S + r)/2 - (?2,^J7^^J}[M] 

-{4:Tt)-"^/^6-^ tj:v{q2,aaLbb -q2.,abLab - 6g2,mm 5" + 6^2,00} [9M] . 

Proof. Define ord(g2.jj) = 0, ord(£;) = 2, ovd{Rijki) = 2, ovd{F) = 2, ord(L) = 1, 
and ord(S') — 1. Increase the order by 1 for each explicit covariant derivative which 
is present. Dimensional analysis then shows the interior integrands in the formula 
for an{Q2, D, B) are homogeneous of order n + 2 while the boundary integrands are 
homogeneous of degree n + I. We use H. Weyl's theorem to write a spanning set 
for the set of invariants and express: 

(4.6) aMQ2,D,B+) = -(47r)-"/2 try(6ig2,^»/2)[M] 

(4.7) a-i(Q2,D,Bt) - -{ATT)-^"'-^^'^^-HTv{ciq2,aa/2 + C2q2,^,n)[dM]. 

(4.8) aQ{Q2.D,B+) = -(4^)-"/26-i{try((72.^(662i? + &3t)/2 

- biq2AjPi])[M] + trv{czq2,aaLbb + CAq2,abLab + C5q2^,nmLaa 
+ CGq2^,nmS + Cjq2^aaS + Csq2,mrn;rn + Cgq2,aa:m)[dM]} . 
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Product formulas then show the constants are independent of the dimension m; 
these invariants form a basis for the integral invariants and are uniquely determined 
for m large. A word of explanation for the formula in equation (4.8) is in order. 
We can integrate by parts to replace the interior integrals q2,ij;ij and q2,ii:jj by 
boundary integrals of (72,mm;m, <?2,am;a, and q2,aa;m- Since q2Mm:a[dM] = 0, we 
use Lemma 4.3 to omit the variable q2.am:a- If we take DM empty, the boundary 
condition plays no role and TZkiik[M] = (sec equation (4.10) below). We use 
Theorem 3.3 to see 6i = 62 = ^3 = ^4 = 1; this completes the proof of assertion (1) 
and the first part of assertion (3). 

The (72, am;* Variables do not appear in equations (4.6), (4.7), and (4.8). Thus we 
may take a variation with 92, am = on dM. This is an essential simplification since 
it means N{q) = g"^"^{Q)^^'^dm- Thus the boundary conditions do not involve any 
tangential derivatives. We set S{g) — g"^"^{g)^'^S. Then the boundary condition is 
preserved; VNie)+S{g) = g"'"'{g)'/\V^ + S). We have 9^5(^)1,^0 = -<Z2,mm5/2. 
By Lemma 2.4, a_i(Q2,£',S+) == -{A7T)-^''^~^')/^A-'^trv{q2.aa/2)[dM]. This shows 
ci — 1 and C2 ~ and completes the proof of assertion (2). 

We use Lemma 2.4 and Theorem 4.1 to see: 

aoiQ2,D,B+) = 
(4.9) - (47r)-("-i)/26-i{try(<72.,(6£; + t)/2 - q2,^,p^J + nkiik)[M] 

+ irv{q2Ma{Lbb + 65*) - 2q2MbLab + "^.Caa - 6<72,mmS') [9M] } 

We have that 

TlkiiklM] — {q2M;ki — q2,ii;kk)[M] = ( — <72, am:a + q2,aa;m)[dM] 

= {Laaq2,rnm — Labq2,ab + <Z2, aa;m) [9M] 
2Caa[dM] ^ (2(72, am;a - <?2,aa;m - (72, mm^aa) [9M] 

= { — iLaaq2,rnm + 2Labq2,ab — q2 ,aa;m)[d M] . 

We use equation (4.10) to compare equations (4.8) and (4.9). This shows 

C3 = 1, C4 = -1, C5 == -2, Cq =: -6, C7 ==6, Cs = 0, Cg = 0. D 

§5 Operators of Dirac Type 

In this section, we study the invariants an{HA,A^), where M is a closed man- 
ifold, A is an operator of Dirac type, and i7 is a smooth endomorphism; we refer 
to Branson and Gilkey [4] for a discussion of the case of manifolds with boundary. 
We begin with a technical result: 

5.1 Lemma. Let A = j'^diy — ip be an operator of Dirac type and let D — A^ be the 
associated operator of Laplace type. Let D — D{g, V, E) and let H G C°° End(y). 

(1) Lo^ = .9.m(-7"5.7" + i^r + r^ + .9"'T,p'0/2. 

(2) Let (/) :— ip + "/'^lo^. Then A = j'^S/u — (j), and (f) is invariantly defined. 

(3) 7»u + 7j;» = 0. 

(4) E = -^^+jt'df,i:-g''^{d^LJ,+LJ,uj^-Lj„r,^^) ^ --fajF,j/2 + j,(f,,,^c^\ 

(5) Let qi,i := -Hji/2, and Qq := -H(p - H.iji/2, then HA = Qi + Qo- 



(4.10) 
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Proof. We compute 

a^ = - 7''a^7^ + 7^-0 + V7^, and b = -V'^ + 7^9^ V'- 

Assertion (1) and the first assertion of (4) follows from Lemma 2.3. We prove 
assertion (2) by computing: 7*Vi — (j) = ^di + Yuii — (j) = 7*9i — ip. We choose a 
system of coordinates and a local frame so that r(xo) = and so that u'(xo)=0. 
Then at xq, we have: 

^-(7'^V,-0)(7^V^-^) 

- (7^7^ + 7^7")/2V,V^ + (7^7^p - H" - 7»V. 

- 7"'/';. + 7 Vf^.M/2 + 0' 
= -5"^V,V^ - E. 

We equate coefficients to derive the second part of assertion (4). Choose a coor- 
dinate system centered at xq so (/^,y = (5^^ + 0(|a:p). We use [3, Lemma 1.2] to 
see that we can choose a local frame for V so 9^7'^ (xo) = 0. Then we have that 
^^(•^o) = ("07^1 + 7ai'0)(2;o)/2. We prove assertion (3) by computing at xq- 

2(7i^;p + 7^;"^) = ['/'7m + 7a'V', li^] + [iplu + lui^, 7^1 
=-07^71^ + 7m "071^ - lui^lt, - lult^i^ 

+ i^lult, + Ivi^lt, - 7m "071^ - llilvTp 

^ipS^^ - S^^ip = 0. 

If we set qi^i = —H^i/2, then Qi = H^iVi + iJ:i7i/2 by Lemma 2.4 since 7^;^ = 0. 
We must therefore take Qq = —H(f> — H.^i'ji/2. D 

The following theorem now follows from Theorem 3.1, Theorem 3.2, and from 
Lemma 5.1. 

5.2 Theorem. Let A = j^di, — ijj be an operator of Dirac type on C'°^'{V) over a 
closed manifold M. Adopt the notation of Lemma 5.L 

(1) ao{HA,A^) - {4TT)-"-/Hrv{Qo}[M]. 

(2) a2{HA,A^) = (4^)-™/26-itry{Oo(T + 6£;) - F,^T,j}[M]. 

(3) ai{HA, A^) == (47r)-'"/2360-i try{go(60£;;fefc + 60tE + 180E^ + Ur-kk 

+5t2 - 2|p|2 + 2|i?|2 + 30FyFy) + 8F,j,kT,j.,k + 8F,j,kqi,kF^, 
~8Fij.kFijqi^k — 4:Fij.jTik:k — 4:Fij.jqi,kFik + 4:Fij.jFikqi^k 

''^^-^ij-^jk^ki I ^ ^ -^ij kn -^ ij ^ kn 1 opjk-t^jn^kn ^^"^^ kn-^kn 

+60E,kqi,kE - 60E,kEqi,k - 30EF,,T,, - 30ET^jF,j}[M]. 
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5.3 Remark. The first two authors [3, Theorem 2.7] computed oq and 02 for / 
scalar and showed that 

(5.4) ao{fA,A^) = - (4^) ""/^ try (/,/.) [M], and 

(5.5) a2{fA,A^)=-{4TT)-"'/%-HTv{f{<pT + 6^E-F,j,j-f,)}[M]. 

For scalar /, Tri2(/.,7,e-*^') = TiL^iAfe-'^'-fAe-'^") == so an{f■^-f^,A^) - 
for all n and we may replace Qo by —fcf) in performing our computations. It then 
follows that Theorem 5.2 (1) agrees with equation (5.4). We see Theorem 5.2 (2) 
agrees with equation (5.5) by using Lemma 2.5 and integrating by parts: 

-tMF^,^^,)[M] = trv((/7,);. - (/7.);,)F,,[M]/2 
= - try(/7jFy;4M]) = try(/7,Fy;,)[M]. 
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